#set up goodness of fit
gm <- modelsummary::gof_map
gm$clean[1] <- "$N$"
gm$clean[7] <- "Adjusted $R^2$"
gm$omit[c(2:6, 8:nrow(gm))] <- TRUE
